Attractive forces between circular polyions 
of the same charge 
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We study two models of ringlike polyions which are two-dimensional versions of simple models 
for colloidal particles (model A) and for rodlike segments of DNA (model B), both in solution 
with counterions. The counterions may condense on Z sites of the polyions, and we suppose the 
number of condensed counterions on each polyion n to be fixed. The exact free energy of a pair 
of polyions is calculated for not too large values of Z, and for both models we find that attractive 
forces appear between the rings even when the condensed counterions do not neutralize the total 
charge of the polyions. This force is due to correlations between the condensed counterions and in 
general becomes smaller as the temperature is increased. For model A divergent force may appear 
as the separation between the rings vanishes, and this makes an analytical study possible for this 
model and for vanishing separation, showing a universal behavior in this limit. Attractive forces 
are found for model A if the valence of the counterions is larger than one. For model B, no such 
divergences are present, and attractive forces are found for a finite range of values of the counterion 
valence, which depends of Z, n, and the temperature. 
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I. INTRODUCTION 

The interaction between colloidal particles in suspen- 
sions are a statistical mechanical problem which has at- 
tracted much attention from both theoretical and exper- 
imental points of view for more than one century now, 
but there are still issues which are not totally understood 
JjJ . Besides the difficulty of properly taking into account 
the long range Coulomb interactions, it is necessary to 
consider the fact that each colloidal particle in solution 
may contribute with a quite large number of counterions, 
and their contribution to the thermodynamic properties 
of the suspension is quite important ||] . One particularly 
important issue related to these systems is the possibil- 
ity of phase separation, and although it has been argued 
that phase separation may occur even in situations where 
repulsive interactions between pairs of colloidal particles 
are present ||], recently it was found that for aqueous 
solutions with monovalent counterions no phase separa- 
tion is expected and this seems to be consistent with 
conclusions drawn from simulations ||. For multivalent 
counterions, however, correlations between counterions 
which are condensed on the colloidal particles may pro- 
duce short range attractive effective interactions between 
the particles, even if their bare charge is not totally neu- 
tralized ||. This interesting phenomenon of effective at- 
tractive interactions between objects with charges of the 
same sign has attracted much attention recently, partic- 
ularly for systems of rodlike polyions, such as segments 
of DNA § 

The theoretical approaches used so far to address the 
problem of attraction between rodlike polyions of the 
same charge may be classified in two groups. The first 
formulates the system as a field-theoretic problem and 



then studies its properties in a high temperature approx- 
imation ||]. The second approach, which is also approx- 
imate, supposes the system of polyions and condensed 
counterions to form a Wigner crystal at zero tempera- 
ture, and then proceeds considering the effect of temper- 
ature on the system |^-[D|. Some time ago, a very simple 
model has been proposed for an isolated pair of rodlikc 
polyions in a solution with counterions which may be ex- 
actly solved, with the aid of a computer, for not too large 
rods [^0) . In this model, the only effect of the condensa- 
tion of a counterion on the site where an ion is located 
is to rcnormalize the local charge. The valence of each 
counterion is a, so that each counterion has a charge aq, 
and the polyion consist of a linear array of Z charges 
—q with a distance b between first neighbors. The so- 
lution has an effective dielectric constant D. The mean 
value of the number n of condensed counterions may be 
determined by the Manning criterion as 
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where 
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is the Manning parameter. For sufficiently low temper- 
ature (£ > 1/a), counterions will therefore condense on 
the polyions. It should be remarked that the expres- 
sion above for the number of condensed counterions is an 
approximation for an isolated polyion and in the limit 
Z — > co. One may wonder if it is still valid for finite 
values of Z and, supposing that the Manning value for 
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n is the mean value of a distribution of numbers of con- 
densed counterions it is relevant to ask if this distribution 
is broad or narrow. It turns out that the expression (|l|) 
is actually a good approximation even for rather small 
values of Z and that in usual physical conditions the dis- 
tribution of values of n is quite narrow |lq| . Now if we 
consider two polyions separated by a distance d, we may 
use expression (JlJ) in the limits d — > oo and d — > 0. The 
total number of counterions condensed on both rods in 
the limit d — > oo will be 



2n 



(3) 



and if we now consider the other limit d — > (where both 
chains merge into one with Z charges — 2q) the number 
of condensed charges will now be given by 
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and therefore the change in the number of condensed 
counterions is given by 
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and for not too high temperatures it may be a reason- 
able approximation to ignore this change, as is done in 
the simple model we described. 

In this paper we study two models for the interaction 
between two ringlike charge structures which are inspired 
by studies of colloidal particles and of rodlike polyions, 
respectively. The models are similar to the one for rod- 
like polyions described above, and we study the ther- 
modynamic properties of these models by solving them 
exactly for not too large particles. In the first model 
(model A), the bare colloidal particles are supposed to be 
formed by Z point charges —q placed regularly on circles, 
and we consider two of these ringlike particles placed in 
a solution with counterions. The number of counterions 
condensed on each particle is equal to n and we proceed 
calculating the force between them as a function of the 
separation between the rings, for fixed values of the va- 
lence a of the counterions and the Manning parameter 
£. The second model is a two-dimensional version for 
a model which was recently proposed for rodlike DNA 
segments, in which the negative charge of the polyion 
is supposed to be concentrated at the axis of a cylinder, 
while the counterions are condensed at the surface of this 
cylinder |fl6| . Thus, in model B we each ringlike polyion 
is formed by a point charge —Zq in the center of the 
circle, while n counterions may condense on Z equally 
spaced sites on the circle. The basic point we address 
in our calculation is the possibility of attractive forces 
between polyions of the same total charge, and for both 
models we find that such forces are indeed possible. 

In the next section we present both models in detail 
and the results we obtained for them, which are cither 
analytical or numerically exact. F ina l comments and dis- 
cussions may be found in section III. 



II. DEFINITION OF THE MODELS AND THEIR 
SOLUTION 



We consider two circular polyions of radius r with an 
even number Z of point charges —q which are placed as 
follows: 

• Model A Uniformly distributed on the circle. 

• Model B In the center of the circle. 

These polyions are in a ionic solution of a-valent counte- 
rions, also supposed to be point charges (charge aq). We 
will admit that n of the counterions are condensed on Z 
sites on the circle. In model A, since the counterions con- 
dense on top of the fixed negative charges, the only effect 
of a condensed counterion is to renormalize the charge on 
the condensation site, which will be equal to q(a — 1). In 
model B, the charge at an condensation site on the circle 
is equal to zero if no counterion is condensed on it and 
aq otherwise. The distance between the centers of the 
polyions is 2r + d. The definitions above are illustrated 
in figure [l]. 





FIG. 1. A possible configuration of two ring polyions with 
Z — 8 and n = 3 condensed counterions. Full small circles 
indicate adsorption sites with condensed counterions, whereas 
empty small circles are empty condensation sites. 

The positions of the condensation sites on the rings will 
be referred by two indices j = 1,2, indicating the ring, 
and i = 0, 1, . . . , Z — 1, indicating the position on the 
ring. For model A, the electrostatic interactions between 
he charges may then be described by the Hamiltonian 
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H A = 



2Dr 4^ , 



(6) 



where the sum is restricted to ^ (i',f) and the 

occupation variables <7i j assume the value if no coun- 
terion is condensed on site (i,j) and 1 otherwise. D is the 
dielectric constant and the ratio of the distance between 
two sites and r is given by 



5(i,j;i',j') = 



(7) 



if j=j', 



[ \JJx + 2 — cos9i — cosOii) 2 + (sinflj — sin^/) 2 otherwise; 
where x — d/r and 
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Model B may be described by the Hamiltonian 
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where the interaction between the central charge and the 
condensed counterions on the same ring was not included, 
since it is constant, and 



S(i) = v / l + (2 + .x) 2 -2(2 + a;)cosl 



(11) 



The partition function Q is calculated summing over 
the possible positions of the condensed counterions 



)a,b = ^2 exp(-/3HA,B), 



(12) 



where the prime denotes the constraint of having n coun- 
terions condensed on each ring. 

It is convenient to define adimensional hamiltonians H 
through 



/3Ha,b = £,'Ha,b, 



(13) 



where £' = [3q 2 /Dr is related to the Manning parameter 
defined in expression (pi) for rodlike polyions. Then 
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and 
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We may define the statistical weights 

Vi = exp (—J ,i = I) 2 ,.. .,Z/2, 



and 
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related to the intra- and interpolyion interactions be- 
tween charges at the condensation sites, respectively, and 



/; = c>q> ( — ) .i = n. l Z/2. 



(18) 



which correspond to the interaction between a central 
charge of a ring and the condensed counterions on the 
other ring. The denominators in the arguments of the 
exponentials are the adimensional distances between the 
interacting charges. For a pair of charges in the same 
poly ion 



(19) 



The expressions for the denominators bi are straightfor- 
ward to obtain from (^|) but will be omitted here for 
brevity, m + 1 = Z 2 /4+Z/2 + 1 is the number of different 
distances between condensation sites located on distinct 
rings. Finally, the adimensional distances between con- 
densation sites and the central charge of the other ring 
are Ci = Si, according to equation (pi]). 

Using these statistical weights, we may rewrite the par- 
tition function as 



Qa.B = y^^A,B{i) 



(20) 
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where the sum runs over the 

Z\ 



n\(Z-n)\) 

allowed configurations of condensed counterions and 
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The exponents UA(ij) and VA(ik) are quadratic polyno- 
mials in a with integer coefficients. For model B we ob- 
tain 



Z/2 m 
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where ubHj) and VB{ik) are integer coefficients multi- 
plied by a 2 and WB{i,l) are negative integer coefficients 
multiplied by Za. The whole set of integer coefficients 
can be calculated with the help of a computer program, 
and then we can obtain the partition function exactly for 
both models and not too large values of Z. 
The force between the rings is given by 



Fa.b = — 



•>A,B 



dd 



(24) 



where 4>a,b = —kgThiQA^B is the Helmholtz free- 
energy. It is convenient to define an adimensional force 
!a,b = (Dr 2 /q 2 )F A ^B, and thus 
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£'Qa,b dx 



(25) 
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For model A, we find that the force between the rings 
becomes attractive for short distances and a > 1. For 
the other cases the force is always repulsive. As x — > 0, 
the force is finite only for a = 1. 




FIG. 2. The force between the rings as a function of the 
distance x for model A. The results presented are for Z = 10, 
n = 2 and f = 2. 

The finite repulsive force observed in the limit x — > 
when a — 1 may be easily understood. In this limit, the 
probability of having no charges condensed on the pair 
of sites (0, 1) and (0, 2) vanishes, since this would lead to 
a divergence in the energy. Thus this pair of sites does 
not contribute at all to the force, and the contributions 
of the other pairs is either positive or zero. This argu- 
ment also explains the decreasing value of the force at 
vanishing distance as n is increased, as may be seen in 
figure [| 

To study behavior of the attractive force for model A 
at short distances (a > 1), we consider the limit of small 
x. In this case, we have 



lim zq = 0, 



(26) 



where zq = exp — £'/ar, is related to the interaction be- 
tween the charges in positions (0,1) and (0,2). All the 
other variables jji and z\ are finite. We can therefore 
define the variable 



v = mm(vA(i, 0)), 



(27) 



where the index refers to the pair of sites (0, 1) and 
(0, 2). We then rewrite the partition function as 



Q A = z^W(l-P), 



(28) 



where 



and 
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(29) 



(30) 



In these equations, we considered the first / configura- 
tions to have an exponent VA(i, 0) = v . P vanishes when 
x — > since va(i, 0) > v for i < I + 1, so we may approx- 
imate 



zlW. 



(31) 

The possible values of va (i, 0), are 1 if there are no con- 
densed counterions in both sites, 1 — a when there is one 
condensed counterion on one of the sites, and (1 — a) 2 
if there is a condensed counterion on each site. As we 
discussed before, the force becomes attractive only for 
a > 1, so we can conclude that the minimum value is 
v = 1 — a. Using the above approximation for the parti- 
tion function we obtain 

1 dW 



x* £W dx 

Taylor expanding the function g(x) = |7^?^T; 9(%) 
go + g\x + ... and keeping only the first term, we get 



(32) 
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(33) 



For the distance xo of vanishing force, we may then find 
the following scaling relation 

xo w ( ) . (34) 
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FIG. 3. The equilibrium distance xo as a function of a — 1 
for £' = 2 (model A). The inset shows data for Z — 8, n — 1 
in the region of a close to 1. The scaling behavior is apparent 
and the numerical fit resulted in a value 0.51 ± 0.01 for the 
slope, in accordance with expression (M). 
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In figure |3| curves of xq as a function of a are shown 
for some values of Z, n, and £, and the collapse of these 
curves at small values of xq may be seen. Also, in the 
inset the limiting behavior found in equation (Q) is evi- 
dent. 

We proceed showing results for model B. In this case, it 
should be stressed that only in the rather unphysical limit 
of infinite temperature we have divergence in the force at 
vanishing values of x. These divergences are present for 
finite temperatures both in model A and in the similar 
model for rodlike polyions [ ^0|Jl7| and are due to pairs 
of sites with nonzero charge with vanishing distance be- 
tween them. The absence of these divergences in model B 
prevents us from performing an asymptotic analysis sim- 
ilar to the one done for model A, which would also prove 
the existence (or not) of attractive forces in model B. 
Simulations performed in a cylindrical version of model 
B found attractive forces at sufficiently low separation 
and enough condensed counterion charges Jig], 
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FIG. 4. The force as a function of the distance for model 
B. The results shown are for Z = 10, n = 2, and a = 1.6. 
Curves for three different values of the inverse temperature 
£' are shown: Dashed line - £' = 2; full line - £' = 1.23; 
dot-dashed line - £' = 0. In the inset, values of the force at 
vanishing x (fo = lim^^o /) are shown as a function of the 
reduced temperature r = l/£'. 

In figure [| the force between the rings for model B 
is shown as a function of the separation x. For suffi- 
ciently low temperatures and not too low or too high 
values of a, attractive forces are found for low separa- 
tion. As mentioned above, for finite temperatures the 
forces are always finite, since in this model the energy of 
the configurations where the pair of sites (0, 1) and (0, 2) 
is occupied diverges as x — > 0. The force at vanishing 
separation fo = lim^^o f{ x ) is shown as a function of 
the reduced temperature r = in the inset of the fig- 
ure, and increases monotonically, so that, as expected, 
it diverges to +oo as the temperature becomes infinite 
(!/£,' — > oo). We obtained similar results for other values 
of Z and n, so that the data in figure ^ are an example 
of a general trend. 
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a-1 

FIG. 5. The equilibrium distance xo as a function of a — 1 
for different values of Z, n, and £' (model B). Full thin line 

- Z = 8, n = 2, and £' = 1; full thick line - Z = 8, n = 2, 
and £' = 0.05; long dashed line - Z = 10, n = 2, and £' = 1; 
short dashed line - Z = 8, n = 3, and f = 1; dot-dashed line 

- Z = 10, n = 3, and f = 1. 

Figure ^| shows results for the equilibrium distance xq 
as a function of the valence a of the counterions for some 
choices of Z, n, and £' . As mentioned above, we were 
not able to perform a simple analytical study of model 
B in the limit of vanishing x, and as expected no sim- 
ple scaling behavior is found for model B in this limit, 
similar to the ones observed for model A and the other 
model proposed for rodlike polyions jl(|[l7]]. Thus, the 
minimum value of a for the existence of an attractive 
interaction varies with Z, n, and and so does the ini- 
tial slope of the curve xq x a. Another feature of these 
curves is the existence of a maximum value of the va- 
lence of the counterions, above which no attractive force 
is found. Thus, as is apparent in figure |5|, for each curve 
there are two values of a for which xq = 0. Results for 
fo as a function of a which will not be presented here for 
brevity also show this feature: The force is attractive in 
the interval between the limiting values of a, and there- 
fore passes through a minimum in some point inside the 
interval. It should be stressed, however, that the larger 
limiting value of a we found for model B assumes val- 
ues which are too large to be physically realizable, but 
since the model itself is rather artificial, one may wonder 
of this effect may not be seen in more realistic models 
such as the one discussed in |l6| ] and in experiments. For 
model A and the model which was proposed for rodlikc 
polymers, no such maximum value of a is found, and the 
curves of xq x a pass through a maximum but approach 
zero asymptotically as a — > oo. 



III. FINAL COMMENTS AND DISCUSSIONS 

In an appropriate range of the parameters, both mod- 
els studied here display attractive interaction forces, 
which may be explained by correlations between the 
counterions condensed on the ringlike polyions. As the 
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temperature is increased, in general the correlations be- 
come weaker and therefore the attractive force decreases. 
In simplified models which were exactly solved before 
0,0, as well as in model A in this work, the attrac- 
tive forces increase without bond as the distance between 
the polyions is lowered, since one or more pairs of sites 
with opposite charges are separated by a distance which 
vanishes. This effect allowed an exact analysis of the 
models at vanishing separation, but also has the draw- 
back of being not very reasonable physically. In partic- 
ular, in the case of model A, it leads to the conclusion 
that it is sufficient to have only one condensed counte- 
rion on each polyion, with a valence slightly larger than 
one, to ensure the appearance of attractive forces. As 
was already found in simulations of cylindrical polyions 
[16 1 and other more realistic models for DNA segments 
[12 1, attraction between like charged polyions is not re- 
stricted to this rather artificial situation where sites with 
opposite charges become very close. In model B we find 
another example of this kind: The forces are finite for fi- 
nite temperature, and a valence which is larger than one 
by a finite amount is needed to ensure attractive forces. 

In both models studied here, when the parameters are 
such that attractive forces are found, the equilibrium dis- 
tance Xq as a function of a has a maximum. At low 
temperature, the maximum is sharp and located close to 
a = Z/n, which corresponds to neutral polyions. As the 
temperature is increased, the maximum becomes broader 
and moves towards larger values of a. With increasing 
temperature, the region of attractive force in the a x x 
plane shrinks and finally disapears as the temperature 
becomes infinite. An example of this behavior is seen 
in the two curves for different temperatures for model B 
with Z = 8 and n = 2 in figure [s| Also, in figure [| one 
notices that for infinite temperature (£' = 0), the force is 
always repulsive. 
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